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Abstract 

In this paper, we study the Cabibbo favored non-leptonic D° decays into K~n + decays. First 
we show that, within the Standard Model, the corresponding CP asymmetry is strongly suppressed 
and out of the experimental range even taking into account the large strong phases coming from 
final state Interactions. We show also that although new physics models with extra sequential 
generation can enhance the CP asymmetry by few orders of magnitude however the resulting CP 
asymmetry is still far from experimental range. The most sensitive New Physics Models to this CP 
asymmetry comes from no-manifest Left-Right models where a CP asymmetry up to 10% can be 
reached and general two Higgs models extension of SM where a CP asymmetry of order 10 -2 can 
be obtained without being in contradiction with the experimental constraints on these models. 
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I. INTRODUCTION 



The Standard Model (SM) has been very successful in predicting and fitting all the 
experimental measurements up-to-date over energies ranging many orders of magnitude l|. 
Unfortunately the SM is only a patchwork where several sectors remain totally unconnected. 
Flavor physics for example involves quark masses, mixings angles and CP violating phases 
appearing in the Cabibbo-Kobayashi-Maskawa (CKM) quarks mixing matrixj^, 3|. These pa- 
rameters unavoidably have to be measured and are independent from parameters present in 
other sectors like Electroweak Symmetry breaking, Quantum Chromodynamics, etc. Other 
sectors remain to be tested like CP violation in the up-quarks sector and even tensions with 
experimental measurements remain to be cleared (see for instance refs.^-7]). 

This is why it is important to find processes where the SM predictions are very well known 
and a simple measurement can show their discrepancy. One of these processes is the rare 
decays and other 'null' tests which correspond to an observable strictly equal to zero within 
SM. So any deviation from zero of these 'null' tests observables is a clear signal of Physics 
beyond SM. This is the case of Cabibbo-Favored (CF) and Double Cabibbo Suppressed 
(DCS) non-leptonic charm decays where the direct CP violation is very suppressed given 
that penguin diagrams are absent j^- 10]. 

Even with the observation of D° oscillation llMl6| and the first signal of CP violation 
in D — > 27r, 2K (Singly Cabibbo Suppressed (SCD) modes) [17H32|. it is not clear that the 



SM 



33 



- |4l| can describe correctly the CP violation in the up quarks sector. It is even more 



difficult as large distance contributions are important and difficult to be evaluated |42|-|46]. 
From the point of view of New Physics (NP), CP violation in CF and DCS modes is an 
excellent opportunity given that it is very suppressed in the SM and it is not easy to find a 
NP model able to produce a reasonable CP violation signal. Thus measuring CP violation 
in these channels is a very clear signal of New Physics. 

Up to now, only D° <B- D° oscillations have been observed and their parameters have 



been measured 
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TABLE I. Direct CP in D non-leptonic decays, from Heavy Flavor Averaging Group HAFG 1 



5l| 



where i ^ or/and y ^ mean oscillations have been observed, while \q/p\ ^ 1 and/or 
7^ are necessary to have CP violation. The theoretical estimations of these parameters [l| 
are not easy as they have large uncertainties given that the c quark is not heavy enough to 
apply Heavy quark effective theory (HQE) (like in B physics) 



471 ] . Similarly it is not light 



enough to use Chiral Perturbation Theory (CP 



cancellations due to the GIM mechanism! 



'h) (like in Kaon physics). Besides there are 
481 ] . Theoretically CP violation in the charm 



sector is smaller than in the B and kaon sectors. This is due to a combination factors: CKM 
matrix elements {\V u bV* b /V us V* s | 2 ~ 1CT 6 ) and the fact that b quark mass is small compared 
to top mass. CP violation in the 6-quark sector is due to the large top quark mass, while in 
the kaon is due to a combination of the charm and top quark. 



Experimental data should be improved within the next years with LHCB 49] and the 
different Charm Factory project 50]. In table the experimentally measured Branching 
ratios and CP asymmetries are given for different non-leptonic D decays. 

In this paper, we study in details the CP asymmetry for the CF D° — > K~n + decay. 



In sect. II, we give the general description of the Effective Hamiltonian describing this 
decay within SM and show how to evaluate the strong phases needed to get CP violating 
observables. These strong phases are generated through Final State Interaction (FSI). In 
sect. Ill, we first evaluate the SM prediction for the CP asymmetry and we show that within 
SM, such CP asymmetry is experimentally out of range. In sect. IV, New Physics models 
are introduced and their contributions to CP asymmetry are evaluated. Finally, we conclude 
in sect. V. 



II. GENERAL DESCRIPTION OF CF NON LEPTONIC D° DECAYS INTO K~ 

AND 7T+ 

In general the Hamiltonian describing D° — > K~n + is given by 



^2 c\ ab sT l c a uTid b + ^ c l 2ab uT l c a sTid b 



(3) 



with i =S, V and T for respectively scalar (S), vectorial (V) and tensorial (T) operators. 
The Latin indexes a, b = L, R and qi, r — (1 ^ 75)9. 

Within the SM, only two operators contribute to the effective hamiltonian for this 



process 



8hll0|. The other operators can only be generated through new physics. 
G F 

H = -j=V* s V ud (ciS7 At c L M7 Al rfi + c 2 w-f^c L s^d L ) + h.c. 
G 



-j=V* s V nd (cA + ^l+Lc. 



(4) 
(5) 



where a\ = c x + c 2 /iV c = 1.2 ± 0.1 and a? = c 2 — C\/N c 

8 
.G F 



-0.5 ±0.1 



color number. For the case D — > K% 

Ao»^K-n+ 

BR 



10l | where N c is the 



4l0| | one has that 



l ^f :gVud 



tdPk 



A\ 



(6) 
(7) 



D 



where BR is the Branching ratio of the process, td is the D lifetime, px is the Kaon 
momentum and mo is the D meson mass. The XJ^> K - and X^° 7r+ can be expressed in the 
following way: 



Xg P3 = if Pl A P2P F P2Ps (m 2 Pl ), A P2Ps = m\ - m\ 



(8) 



where fo and fx are the decay constants for D and K mesons respectively and F® and 
Fq* are the corresponding form factors. These amplitudes have been computed within 
the so called naive factorization approximation (NFA) without including the Final State 
Interaction (FSI). In NFA, no strong CP conserving phases are obtained (and therefore no 
CPV is predicted) but it is well known that FSI effects are very important in these channels 



52 



561 ] . In principle you have many FSI contributions: resonances, other intermediate 



states, rescattering, and so on. Resonances are specially important in this region given that 
;hey are abundant. They can be included and seems to produce appropriate strong phases 



56(. However the other contributions mentioned above have to be included too, rendering 



the theoretical prediction cumbersome. A more practical approach, although less predictive, 



is obtained by fitting the experimental data 



52, 



561 ] . This is the so called quark diagram 



approach. Within this approach, the amplitude is decomposed into parts corresponding to 
generic quark diagrams. The main contributions are the tree level quark contribution (T), 
exchange quark diagrams (E), color-suppressed quarks diagrams (C). Their results can be 
summarized in the following way, for the process under consideration 56| : 

Ad°-+k-it+ = V* s V ud {T + E) (9) 

with 

T = (3.14 ±0.06) • 10~ 6 GeV 

E = 1.53^-10- 6 -e( 122±2 )° l GeV (10) 
where in NFA they can be approximately written as 

T ~ ^a.Uml - ml)F DK (ml) (11) 

E ~ -^=a 2 f D (ml - ml)Ft{m 2 D ) (12) 

In the rest of this work we are going to use the values obtained by the experimental fit, 
given in eq. (TlOl . 



III. CP ASYMMETRY IN D° -»• K~ir + WITHIN SM 



In the case of CF (and DCF) processes the corrections are very small (see diagrams in 



fig.l and fig. 2) and are generated through box and di-penguin diagrams 
section, we shall evaluate these contributions. 



57 



59|. In this 
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FIG. 1. Feynman diagram for CF processes: Box contribution. 



The box contribution is given as 



59 
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(13) 
(14) 



where 



b x = X^ d f( Xu , x D ) (15) 

= v; d v ud (v: s v ud f ud + v; s v cd f cd + v t * s v td f td ) 
+v; s v us (v: s v ud f us + v: s v cd f cs + v*v td f ts ) + v; b v ub (v* s v ud f ub + v c * s v cd f cb + v*v td f tb ) 

= Vastus Wcs^cd (fcs ~ fed — fus + fud) + ^tl^td (fts ~ ftd ~ fus + fud)} 

+v; b v ub [v; s v cd (u - f cd - f ub + f ud ) + v*v td (f tb - f td - u + f ud )} (16) 



with A^ D , = VirpVnn', = Vj} D Vu> D) U = u, c, t and D = d, s, 6, x q = (m q /m w ) 2 and 



fuD = f{xu,X D ) 



61] 
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Numerically, one obtains 



b T ~ 3.6 • HTV 



7„0.07-i 



(17) 



The quark masses are taken their values at m c sea 



to the Lagrangian is the dipenguin and it gives 
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62] 



given in [![. The other contribution 
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FIG. 2. Feynman diagram for CF processes:di-penguins contribution. 



AW = [\°Eo(*d)] K d E (xu)] s lfl T a d L {jsTU - u lu T a c L 

= -^fP9sl,T a d L {sTU - d»d v ) u lu T a c L (18) 

I67H 

p g = [Xf u E (x D )] [\V d E ( Xu )] = [V; S V US (E (x s ) - E (x d )) + V; b V ub (E (x b ) - E (x d ))} 
[V cd V c * s (E (x c ) - E (x u )) + V td V* (E (x t ) - E (x u ))} (19) 



where T a are the generator of SU(3)c- Numerically, p g ~ —1.62 • e °- 002i and the Inami 
functions are given by 

E (x) = — — — - \x(l - x)(18 -llx-x 2 )- 2(4 - 16a: + 9x 2 ) \og(x)] (20) 
12(1 — x) 4 

The operator O can be reduced as 

O = s^T a d L (g^D - 8»d v ) u lv T a c L = s^T a d L D (uYT a c L ) + s$ T a d L u$ T a c L 
= -q 2 s lfl T a d L u^T a c L - (m s sT a d S -P + m d sT a d s+P ) ■ (m c uT a c s+P + m u uT a c S - P ) 
-q 2 s^T a d L uYT a c L - m s m c sT a d L uT a c R - m d m u sT a d R uT a c L 
—m s m u sT a dLuT a cl — m d m c sT a dRuT a CR (21) 

where q 2 is the gluon momentum and N is the colour number. This expression can be 
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simplified using the fact that 



s^T a d L uYT a c L = \{0i- ^0 2 

sT a d L uT a c R = -^,iCRvrfd L - ^sd L uc R 
sT a d R uT a c L = ~s-f^c L u^d R - -^sd R uc L 



sT"d L nT"c L --sc L ud L - Y^sa^CLua^di - —sd L u<- L 



1_ _ 1 _ _ 1 
sT a d R uT a c R = --sc R ud R - —sa^ u c R ua^d R - —sd R uc R (22) 



Once taking the expectation values, one obtains 



(O) = -q 2 (s^T a d L u^T a c L ) - m s m c (sT a d L uT a c R ) - m d m u (sT a d R uT a c L ) 
-m s m u (sT a d L uT a c L ) - m d m c (sT a d R uT a c R ) 

~~2\ N*J Xd ° k ~ + — { 1 ~n) Xd ° k ~ + 8N^ s mDXK ^ + (23) 



Hence, one gets for the Wilson coefficients 



G F m 2 w G F a s r :1 



Aai = r- — ,A _ 



2.8 • IQ-V - 0044 



q ( 1 \ m c m s f 1 



2 V 1 iV 2 / 4 l 1 



Gpmly G F a s 5m d m 2 D 
Aa 2 = = b x ■= —Pa 

V2 irW*V ud A^V CS V: S 8Nm s 

~ -2.0 • l(T 9 e acm (24) 

where to obtain the last result it has been used the fact that for the decay D° — > K~tt + , 
one can approximate q 2 = (p c TPu) 2 = (p s ±Pd) 2 — (pd ~ Pn/2) 2 = ijn 2 D + m 2 K )/2 + 3m 2 r /4, 
by assuming that p c ~ po and p u ~ p,r/2 and as — 0.3. It should be noticed that the box 
contribution is dominated by the heavy quarks while the penguin is by the light ones. The 
direct CP asymmetry is then 

\A\ 2 - \A\ 2 2|r|sin(0 2 -0 1 )sin(o^))) 10 

Acp ~ \a\ 2 + \a\ 2 = \TT^f " L4 ' 10 (25) 

with r = E/T, — >• a« + Aaj = + |Aaj| exp[iA0j] and fa ~ Aa^sin A0j/aj and ag is the 
conserving phase which appears in eq. ffTOj) . 
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IV. NEW PHYSICS 



With New Physics, the general Hamiltonian is not only given by C?i.2- The expressions 
of the expectation values of these operators can be found in the appendix. It is important 
to notice that as expected only two form factors appear, namely Xk-tt+ anc ^ Xd°k-- This 
is important to take into account the FSI interactions as the first one is identified as E 
contribution and the second one is identified as T contribution. In the next subsections, we 
shall calculate the Wilson coefficient for different models of New Physics. The first case will 
be assuming to have extra SM fermion family. The second example will be to compute the 
CP asymmetry generated by a new charged gauge boson as it appears for instance in models 
based on gauge group SU{2) L x SU{2) R x U(1)b-l an d our last subsection is dedicated to 
the effects CP asymmetry coming from new charged Higgs-like scalar fields, applying to two 
Higgs extension of the SM (type II and type III). 



A. Contributions to Aqp from extra SM fermion family 



A simple extension of the SM is the introduction of a new sequential generation of quarks 



and leptons (SM4). A fourth generation is not exclude by precision data[63l-7o||. Recent 



reviews on consequences of a fourth generation can be found in 

The B — > Kit CP asymmetries puzzles is easily solved by a fourth generation 
with a mass within the following range |82]: 

400 GeV <m UA < 600 GeV. 



82 



24 1 



(26) 



The value of SM4 parameters compatibles with the high precision LEP measurements 
66l . |69J are 

1 



m 



11,4 



1 + 5 ln Tl5GeV 



x 50 GeV 



V ud X\V UA d\ < 0.04 



(27) 
(28) 



where V is the CKM quark mixing matrix which is now a 4 x 4 unitary matrix. The direct 
search limits from LEPII and CDF 85H871] are given by: 



m U4 > 311 GeV 



(29) 



m d4 > 338 GeV. 



Direct search by Atlas and CMS coll. have excluded m^ 4 < 480 GeV and m q4 < 350 GeV 
88-1901] . above the tree level unitarity limit, m U4 < v ~ 504 GeV. But SM4 is far 



to be completely understood. Most of the experimental constraints are model-dependent. 
For instance it has been shown in 9Jj that the bound on m U4 should be relaxed up to 



m 



"4 



> 350GeV if the decay u<± — > ht dominates. The recent LHC results which observe 



an excess in the H — y 77 corresponding to a Higgs mass around 125 GeV 92|, |93| seems 



to exclude the SM4 scenario 



94| but this results is based on the fact that once we include 



the next-to leading order electroweak corrections, the rate cr(gg —y H) x Br(H —y 77) 
is suppressed by more than 50% compared to the rate including only the leading order 



corrections 



94ll99j . This could be a signal of a non-perturbative regime which in SM4 



can be easily reached at this scale due to the fourth generation strong Yukawa couplings. 
Therefore, direct and model-independent searches for fourth generation families at collider 
physics are still necessary to completely exclude the SM4 scenario. 

The CP asymmetry in model with a fourth family is easy to compute as the contributions 
come from the same diagrams in the SM with just adding an extra U4 = t' and d± = V . 
Similarly in ref. 90], it has been found that new CKM matrix elements can be obtained (all 
consistent with zero and for my = 600 GeV) to be 



Su 

\V« d \ 
tan 

\Vcbl 



\V u y \ = 0.017(14), S 24 



I V c y 



0.0084(62) 



w, 



lb' 



0.07(8) 



\V t ' 



t's\ 



v.. 



V, 



ad 



0.01(1), \v fb 



C14 C14 C14C24 C14C24 

0.07(8), \V V y\ = 0.998(6), \V tb \ > 0.98 



Sl3 



C14 



5 13 = 7 = 68° 



IC13C24S23 - U 13 UuU 24 \ ~ C13C24S23 



(30) 



The two remaining phases (<f)u and ^24) are unbounded. Thus the absolute values of the 
CKM elements for the three families remain almost unchanged but not their phases. From 
these values one obtains 

S13 = 0.00415, su = 0.225, s 23 = 0.04, s u = 0.016, s 24 = 0.006, s 34 = 0.04 (31) 

For a 4th sequential family the maxima value for the CP violation is obtained as 

A CP ~ -1.1 • 10" 7 

(32) 



where one uses \V u y\ = 0.06, \V r y \ = 0.03, \Vty\ = 0.25 



'14 



-2.9, 



'24 



1.3 
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This maximal value is obtained when the parameters mentioned above are varied in a 
the range allowed by the experiential constrains, according to eq. [30] in a 'three sigma' 
range. The phases are varied in the whole range, from — 7r to n. Thus one can obtain 
an enhancement of thousand that may be large but still very far from the experimental 
possibilities. 



B. A new charged gauge boson as Left Right models 



In this section, we shall look to see what could be the effect on the CP asymmetry 
coming from a new charged gauge boson coupled to quarks and leptons. As an example of 
such models, we apply our formalism to a well known extension of the St andard Model based 



on extending the SM gauge group including a gauge SU(2) R 104j | . So now, our gauge 
group defining the electroweak interaction is given by SU (2) L x SU (2) R xU (l)s_r^This_SM 



CMS 



115 



extension has been extensively studied in previous works (see for insta nce refs. Il05l - ll09| ) 
and t heir param eters have bee n strongly constrained by experiments (lL IllQ - 1 14 1 . Recently, 



116 and ATLAS 

rT 

the Wr gauge boson mass 119 . 



117 



118 



at LHC have improved the bound on the scale of 
The new diagrams contributing to D — > Kit are similar to 
the SM tree-level diagrams with Wl is replaced by a Wr. These diagrams contribute to the 
effective Hamiltonian in the following way assuming no mixing between Wl and Wr gauge 
bosons : 



H 



LR 



( 



9r™w 



\Pi \9Lm WR 
Gp 

71 



VRc S V Rud (c^sj^CRUj^dR + c^wy^CRS^dn) + H.C. 
V*V vd (c 1 1 + C20 2 )+li.C. (33) 



where gi and g R are the gauge SU (2) L and SU(2) R couplings respectively. m w and m WR are 
the SU{2) L and SU{2) R charged gauge boson masses respectively. Vr is the quark mixing 
matrix which appears in the right sector of the lagrangian similar to the CKM quark mixing 
matrix. This new contribution can enhance the SM prediction for the CP asy mme try but 
still it is suppressed due to the limit on Mw R which has to be of order 2.3 TeV jll9j | in case 
of no-mixi ng Left r ight models. 



In refs. 



125| it has been shown that the mixing between the left and the right gauge 



bosons can strongly enhance any CP violation in the Charm and muon sector. This LR 



11 



mixing is restricted by deviation to non-unitarity of the CKM quark mixing matrix. The 
results were that the Left-Right (LR) mi xing angle called £ has to be smaller than 0.005[121| 
and right scale M R bigger than 2.5 TeV 1191 ]. If the Left-Right is not manifest (essentially 
that g R could be different from at Unification scale), the limit on Mr scale is much less 
restrictive and the right gauge bosons could be as light as 0.3 TeV [1221 ]. In such a case, 
£ can be as large as 0.02 if large CP vi olation ph ases in the right sector are present 107 ] 



still compatible with experimental data 



123 



- Il25| . Rece ntly, preci sion measurement of the 



126 



127] put model independent 



muon decay parameters done by TWIST collaboration 
limit on £ to be smaller than 0.03 (taking = g R ). Let's now compute the effect of the LR 
mixing gau ge boson on our CP asymmetry. So first, one defines the charged current mixing 



matrix 



120 ] 





H 


\ W R ) 





cos £ — sin £ 

e tul sin £ e JW cos £ 





(34) 



where W\ and W% are the mass eigenstates and £ ~ 10 2 . Thus the charged currents 
interaction part become 

£ ~ ~^f ltl {qlVPl + g R Z,V R P R ) DW{ - -±=U^ (-g L £VP L + g R V R P R ) DW} (35) 

where V = Vqkm an d V R = e luJ V R . Once one integrates out the W\ in the usual way and 
neglecting the W2 contributions given its mass is much higher, one obtains the effective 
hamiltonian responsible of our process: 

IGV L s 7 „ ( V*P L + °*tV R *P R ) c uY (VP L + —HVP R ) d 
L V 9l J cs V 9l J ud 

+ h. c. 



n 



cff. 



V2 



9l 



< 2 ( V*P L + ^V R *P R ) cpufff ( VP L + y^yPn ) <L 



ad 



(36) 



where a, (3 are color indices. It is easy to check that taking the limit £ — > 0, one obtains 
eq.flH]) with the only difference comes from the C2 terms, the Fierz transformation has been 
applied. The terms of the effective Hamiltonian proportional to £ are: 

AH eS ~ [ Cl s lt y*c L uYV^d R + Cl s^V R *c R uYV ud d L 

V 2 9l 

C2S a l ^V* s c L pu^V R d d Ra + c 2 s a j IJi V R *c R pUpj' J, V u dd La \ + h. c. (37) 
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The contribution to the amplitude proportional to £ is then given by: 
iG F g R 



AA 



V2 Ql 



2 Y D ° 



-c2v; s v u R d (2 X D °x^ + + ±x$ K . 

iG F g R 



■ad 



+ c l V c fV l 

- c 2 v c fv ud (2 X D °x^: w+ + -x D/0K 



2 Y D ° 



1 

N' 



/- k (V*V u R d - V c fV ud ) (axX£ K - + 2 X D °a 2 X^ 
V 2 9l v 

= - 9 -i {v c fv ud - v; s v u R d ) (t - 2 X D °E 

wherex ?r+ and are defined as 



(38) 



X" 



ml 



(m c - m s )(m u + m d ) 



x D ° 

The CP asymmetry becomes 



??? 



(m c + m u )(m s - m d ) 



(39) 



.4 



CP 



1 + 2 X D ° ) Im {V c fV ud - V*V*) Im(r) 



V*V ud \l + r\ 2 
with r = E/T. For a value as large as £ 



(40) 



10~ 2 the asymmetry can be as large as 0.1. Also, 
we should notice that to obtain this results, it has been used the fact that the chiralities 
don't mix under strong interactions, if the quark masses are not taken into account. This 
is approximately the case in the evolution of the Wilson coefficients from my/ to m c as the 
quark in the loop are the down quarks contrarily to process like 6—7-57 where the quarks 
in the QCD correc tions are the up quarks and in that case, a strong effect from top quarks 



could be expected 



1 2 8N 1 3 lj ] . In our case, as a first approximation, the QCD corrections to 



the Wilson coefficient coming from the running of the renormalization group from mw to 
m c can be safely neglected. 



C. Models with Charged Higgs contributions 



Our last example of new physics is considering contribution to the effective Hamiltonian 
responsible of the D° — > K~tt + process due to a new charged Higgs fields. The simple 
SM exte nsions wh ich include new charged Higgs fields are the two Higgs doublet models 
(2HDM) 132l . [133(1 . Usuall y, it is used to classify these 2HDM in three types: type I, II or 



III (for a review see ref. 1341 ]). In 2HDM type II models (like Minimal Supersymmetric 
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Standard Model), one Higgs couples to the down quarks and charged leptons and the other 
Higgs couples to up type quarks. LEP has performed a Direct search for a charged Higgs 
in type II 2HDM and they obtained a bound of 78.6 GeV 135] . Recent results on B — ^ tv 
obtained by BELLE [5] and BABAR [(J h ave strongly improved the indirect constraints on 
the charged Higgs mass in type II 2HDM 1361 ]: 



m H + > 2A0GeV at 95%CL 



(41) 



2HDM type III is a general model where both Higgs couples to up and down quarks. Of 
course, this means that 2HDM type III can induce Flavor violation in Neutral Current and 
thus it can be used to strongly constrain the new parameters in the model. We shall focus 
our interest to the two Higgs doublet of type III as the other two can be obtained from type 
II taking some limits. In the 2HDM of type III, the Yukawa Lagrangian can be written as 



1371 . 



138) : 



- Q a fL [Yfc ab H h : + e u fi H°i\ u lR + H.c. 



(42) 



where e & is the totally antisymmetric tensor, and e^- parametrizes the non-holomorphic 
corrections which couple up (down) quarks to the down (up) type Higgs doublet. After elec- 
troweak symmetry breaking, £y gives rise to the following charged Higss-quarks interaction 
Lagrangian: 

r e ff —t, r H±LRcS p a _L,-7 r H±RLcS P a (ai\ 
L H± - UfL Ufd . F R di + UfL Ufd . F L di , (43) 



with [138) 



tH ± LReft 
Ufdi 



-iif± RL eff 
urdi 



sm P V fi C^Sji - 4 tan fi\ , 

3 

cos/3 



m 



(44) 



Here v u and are the vacuum expectations values of the neutral component of the Higgs 
doublets, V is the CKM matrix and tan(3 = v u j Vd- Using the Feynman-rule given in Eq. fj4*3|) 
we can compute the effective Hamiltonian resulting from the tree level exchanging charged 
Higgs diagram that governs the process under consideration namely, 

n eff = ^v: s v ud c?0jl)q*m, (45) 



1=1 
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where Cf are the Wilson coefficients obtained by perturbative QCD running from M#± 
scale to the scale \x relevant for hadronic decay and Qf are the relevant local operators at 
low energy scale // ~ m c . The operators can be written as 

Qf = (sP R c)(uP L d), 
Qf = {sP L c){uP R d), 
Qf = (sP L c)(uP L d), 

Qf = (sP R c)(uP R d), (46) 
And the Wilson coefficients Cf , at the electroweak scale, are given by 



c; 



C " = WOT (?-' * (**.- 4 *» ) ( E - * V V (=«* - 4 tan / j) ) 



fc=l ^ ' ' X J=1 

We now discuss the experimental constraints on the e?- where q = d,u. The flavor-changing 
elements e^- for i ^ j are strongly constrained from FCNC processes in the down sector 
because of tree- level neutral Higgs exchange. Thus, we are left with only ej 1; ef 2 - Concerning 
the elements & we see that only e^^e^ can significantly effects the Wilson coefficients 
without any CKM suppression. Other efj terms will be so small as the CKM suppression will 
be of orders A or A 2 or higher and so we neglect them in our analysis. One of the important 
constraints that on e|- where q = d,u can be obtained by applying the naturalness criterion 
of 't Hooft to the quark masses. According to the naturalness criterion of 't Hooft, the 
smallness of a q uant ity is only natural if a symmetry is gained in the limit in which this 
quantity is zero 138| . Thus it is unnatural to have large accidental cancellations without a 
symmetry forcing these cancellations. Applying the nat uralness criterion of 't Hooft to the 
quark masses in the 2HDM of type III we find that 138 ] 

K(d)<4 | < I'VijI max [m d ^ Ui ),m d ^ Uj )] . (48) 
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(47) 
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FIG. 3. Constraints on e 22 . Left plot corresponding to tan/3 = 10 while right plot corresponding 
to tan/3 = 100. 



which leads to 



< 





max [ 






V u (d) 





(49) 



Clearly from the previous equation that e" 1; ef 1: e 22 will be severely constrained by their small 
masses while e 22 will be less constrained. Clearly from Eq. (l49j) . the constraints imposed on 
are tan/3 dependent. We now apply the constraints imposed on the real and imaginary 
parts of €32 corresponding to two different values of tan/3 namely for two cases tan/3 = 10 
and tan/3 = 100 using Eq. (l4"9"j) . In Fig. (j3J) we show the allowed regions for the two cases. 
Clearly the constraints are sensitive to the value of tan (3 where the constraints are weak for 
large values of tan (3. Since C± and Cf are proportional to e 22 thus they will be several order 
of magnitudes larger than C 2 and . In fact this conclusion can be seen from Eq. (j47p 
and thus in our analysis we drop C 2 and C% ■ Now possible other constraints on e 22 can 
be obtained from D — D mixing, K — K mixing. For K — K mixing, the new contribution 
from charged Higgs mediation corresponding to top quark running in the loop will be much 
dominant than the contribution in the case where the charm quark runing in the loop. This 
is due to the dependency of the contribution on the ratio of the quark mass running in the 
loop to the charged Higgs mass. Thus the expected constraints fr om K — K mixing might 

the constraints on 



be relevant on e| 2 and not on e 22 . In fact, as mentioned in ref. 



138 



eg 2 and are even weak and e^and can be sizeable. By a similar argument we can 
neither use the process b — > S7 nor the Electric dipole moment (EDM) to constraint e 22 . 
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Regarding D — D mixing one expects a similar situation like that in K — K about the 
dominance of top quark contribution. However due to the CKM suppression factors the top 
quark contribution will be smaller than the charm contribution. 



1 . D — D mixing constraints 

We take into accounts only box diagram that contribute to D — D mixing mediated 
by exchanging strange quark and charged Higgs. Other contributions from box diagram 
mediated by down or bottom quarks and charged Higgs are suppressed by the CKM factors. 
Since SM contribution to D — D mixing is very small we neglect its contribution and neglect 
its interference with charged Higgs mediation contribution. Thus effective Hamiltonian for 
this case can be written as: 




i=l 



(50) 



where Cj, C t are the Wilson coefficients obtained by perturbative QCD running from M H 
scale to the scale \x relevant for hadronic decay and Qi, Qi are the relevant local operators 
at low energy scale 



Q 1 = (uYP L c)(u^P L c), 
Q2 = (uP L c)(uP L c), 
Q 3 = (uYP L c)(u ltl P R c), 
Q 4 = (uP L c)(uP R c), 



(51) 
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where we drop color indices and the operators Qi can be obtained from Qi by changing the 
chirality L R. The Wilson coefficients Cj, are given by 

3 / \\2/3 / \ \ 2 



Ci 



Co 



c 3 



647T 5 



167i 2 m 2 LJ 



sm /? ^ (^,2 - e? 2 tan ^ ( £ sin /3 (^ fc2 - 2 tan ) 
£ sin Z 3 ^ f — ^2 - e% tan /?) ) ( £ cos /3 (—S kl - e£ tan /3 

7=1 yy \ k =i ^ Vu 



64vr 2 



£ sin /3 V 2 * ^5 i2 - tan ^ sin /3 V u (^<5 fc2 - e d k2 tan ) 

7 = 1 fe=l 



a, 



^ cos p V l2 {^5 a - ejf tan ^ ^ cos ^ (^na - e£ tan /3 

£ sin V* (^5 j2 - e% tan p\ V £ sin /3 Vi* ( 

7=1 ^ Ud ' ' \ k=i ^ 



m 2 s I 2 (x s ) ' 
167r 2 m 2 H± 



v d 



$k2 - £ d k 2 tan ^ 



£ cos /3 \/ /2 ^(Jn - ejf tan ^ (j£ cos /9 (^2 - e£ tan ^ ) . 



(52) 



where x s = m 2 /m H± and the integrals are defined as follows: 

x s + 1 -2x s ln(x s ) 

= ~ + 



x, - 1 



fx, - 1) 3 



2 + (x s + l)ln(x s 



( Xs _l) 2 (x s -l) 3 
The Wilson coefficients Cj are given by 



(53) 



C 2 



647T 5 

m 2 I 2 (x s 



^cos/J^^^-e^tan/?^ ^ cos /3 y fe * 2 - e& tan ^ , 

£ cos /3 V£ ^5,- 2 - e£ tan ^ ( sin ^ Vlk (^ fc2 ~ e " 2 tan ^) ) 



lQ7r 2 m 2 1± 



C3 — C 3 , 

C4 = C4. 



(54) 



Our set of operators Qi , Q? and Q 4 given in Eq. ([5T]) are equivalent to their corresponding 
operators given in Refs. [1391 . Il40| while the operators Qi and Q 2 are equivalent to Qq and 
Qi given in th e same re ferences respectively Moreover Q 3 , given in Eq. (l5ip . can be related 



to Q 5 in Refs. 



139, 



140l | by Fierz identity For the rest of the operators, Q 3 and Q4, they 
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are equivalent to Q$ and Q4 in Refs. [1391 . Il40| since their matrix eleme nts are eq ual. Thus 



our Wilson coefficients can be subjected to the constraints given in Ref. 
find that 

r fn H ± 



139 



140| and so we 



Cx I < 5.7 x 1(T 7 
C 2 I < 1.6 x 1(T 7 
C 3 I < 3.2 x 10~ 7 
C A I < 5.6 x 1(T 8 



ITeV 
ITeV 

ITeV 

m H ± 
ITeV 



(55) 



the constraints on C\ — C4 are similar to those C\ — C 4 . As can be seen from Eq.( l55]) the 
constraints on the Wilson coefficients will be strong for small charged Higgs masses. We 
can proceed now to derive the constraints on e 22 using the upper bound on C 2 for instance. 
Keeping terms corresponding to first order in A where A is the CKM parameter we find that, 
for rriH± = 300 GeV and tan /3 = 55 

C 2 x 10 12 ~ 3 ( - 53.6 e d 12 - 12.7 e d 22 + 0.007 



- 12.4 e\% - 53.4 e u 22 * + 0.007 (56) 



While for m H ± = 300 GeV and tan (3 = 500 we find 



Co x 10 



ii 



3.6 



487.1 e d 12 - 115.0 4, + 0.06 



112.5 6^ - 486.7 t u 22 * + 0.007 (57) 



In both Eqs.( l57f56|) we can drop terms proportional to e" 2 * to a good approximation as 
they have small coefficients in comparison to e 22 and also since e^- d with i 7^ j are always 
smaller than the diagonal elements e\f '. On the other hand we know that ef 2 can not be 
large to not allow flavor changing neutral currents and so we can drop terms proportional 
to ef 2 in Eqs.( l57l56l) to a good approximation also, thus we are left with e 22 and e 22 in both 
Eqs.( l57f56|) . Comparing their coefficients shows that e 22 has a large coefficient and thus we 
can drop e 22 terms. An alternative way is to assume that e 22 terms are the dominant ones 
in comparison to the other e*- d terms and proceed to set upper bounds on e 22 . In fact even 
if we consider other Wilson coefficients rather than Co this conclusion will not be altered. 



Under the assumption ef 2 = e 



22 



: 12 



and using the upper bound corresponding to 



m H ± = 300 GeV on C 2 , using Eq.( )55|) . which reads in this case 

I C 2 I < 1.4 x 10" 8 



(58) 
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Clearly from Eqs.f l56|l57f58|) the bounds that can be obtained on will be so loose and 
thus D — D mixing can not lead to a strong constraints on e^. 



2. D q — > tv constraints 



The decay modes D q — > tv where q = d or q = s can be generated in the SM at tree level 
via W boson mediation. Within the 2HDM of type III under consideration, the charged 
Higgs can mediate these decay mode s at tree level also and hence the total branching ratios, 
following a similar notations in Ref. 138] , can be expressed as 



B(D1 



r 2 It/ i 2 

F \ v cq\ 2 r2 

-m T f D m Dq 



57T 



1 + 



m 



Da 



711 



cq * 
R 



'm c + m q ) m r 



s^cq* 



Where we have used 141 ] 



mrc\D q ) 



fD q m 2 Da 



(59) 



(60) 



(m c + m q ) 

Where the SM Wilson coefficient is given by = AGf V^/V?* and the Wilson coefficients 
C c £ and at the matching scale are given by 



ricq 
°i?(L) 



M 



J_ V LR{RL),H± r^ tan/3; 



(61) 



with the vacuum expectation value v ~ 174GeV and r^(- Ri )>- ff can be read from Eq. 
Setting the charged Higgs contribution to zero and fo s = 248 ± 2.5 MeV jl42 |. we find 
that B SM (D+ ->• t+v) ~ 9.5 x 10~ 4 and B SM (Dt ->■ t+v) = (5.11 ± 0.11) x 10~ 2 which 



is in close agreement with the results in Ref. 143l - ll45| . The exp erimental values of these 
Branching ratios are given by B(D^ — > t + v) < 2.1 x 10~ 3 146] while B(D+ — > t + v) = 
(5.38 ± 0.32) x 10~ 2 142] . Keeping the terms that are proportional to the dominant CKM 
elements we find for q = d 



^H ± RL eff 
cd 

-iif± LRcS 
cd 



COS /3 Vxi (—612* tan 0) 

m d 



sin/31/] 



ii 



Vd 



- e\ x tan (3 



(62) 
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While for q = s we find 



rf ±Meff = cos^y 22 (^-e^tan/3 



r£ ± Meff = sin p V 22 ( ^ - <& tan /3 



I'd 



(63) 



Clearly from the last two equations, we need to consider the decay mode — >• t + v to 
constrain e^. For tan/3 = 10 we find that 



r H±^eff xl0 -3 
r ^LHe ff><10 -3 



0.71 - 968.6 e^ 2 



5.3 - 9686.0 e 
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(64) 



Clearly the coefficient of e 22 is one order of magnitude larger than e 22 and for larger tan /3 
one expects to be larger than. However, e 22 is severely constraint by naturalness criterion 
and thus we expect that the term proportional to e^ 2 to be larger and thus in our analysis 
we can drop e 22 term and proceed to obtain the required constraints. We show in Figs.( )4f5l) 
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FIG. 4. Constraints on e 22 from B(D^~ — > T^"f). Left plot corresponding to tan/3 — 200 while right 
plot corresponding to tan /3 = 500. In both cases we take m#± = 200 GeV. 

the allowed regions for the real and imaginary parts of e 22 corresponding to two different 
values of the charged Higgs mass namely, m#± = 200 and rriH± = 300 and for different 
values of tan/3. Our objective here is to show the dependency of the constraints on m#± 
and tan/3. We see from the Figures that, for tan/3 = 500, the constraints become loose 
with the increasing of m#± . This is expected as Wilson coefficients of the charged Higgs are 
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FIG. 5. Constraints on from B(D+ —> t + u). Left plot corresponding to tan/3 = 350 while right 
plot corresponding to tan /3 = 500. In both cases we take m#± = 300 GeV. 

inversely proportional to the square of m#± and thus their contributions to B(D+ — > t + u) 
become small for large m#± which in turn make the constraints obtained are loose. Another 
remark from the figure is that the constraints become strong with the increasing of the value 
of tan/3 which is expected also from Eq. (l6T|) . This in contrast to the constraints derived by 
applying the naturalness criterion where we showed that the constraints become loose with 
the increasing of the value of tan (3. 

3. CP violation in Charged Higgs 



The total amplitude including SM and charged Higgs contribution can be written as 



A = [c! M ^c! M + x *\c»-c»))x^ 

(65) 

with Xp^ Pa = if Pl Ap 2 p 3 FQ 2P3 (m Pi ), ^p 2 p 3 = rn% 2 - mp 3 and x n+ and x D ° are previously 
defined as 

,2 



r<SMj_ 1 r>SMj_ 1 (r<H ^D°r<H\ \ v d° 
c 2 +]^ G i + ^/v^°i ~ x L 'i)) X K-^ 



mi 



X 



(m c - m s )(m u + m d ) 



X 



D° 



III 



D 



(66) 



(m c + m u ){m s - m d ) 

The form of the amplitude, A, shows how charged Higgs contribution can affect only the 
short physics (Wilson coefficients) without any new effect on the long range physics (hadronic 
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parameters) . Thus strong phase will not be affected by including charged Higgs contributions 
while the weak phase will be affected. We can rewrite Eq. (l65j) in terms of the amplitudes T 
and E introduced before in the case of the SM as follows: 



A = V* s V ud (T SM+H + E SM+H ) (67) 



where 



rpSM+H = x 10 -6 „ <fF a SM+H f ^ m 2 D _ ^pDK ^ 

V2 

E SM+H = L5 g x 10 -6 e m°< „ GfgSM+H _ ^pK,^ (6g) 

v2 



where 



,s u + U I C SM + }_ C SM + y , {( .U _ ( .U 



,SM+H 



ax + Acn + x* (Cf -Cf) (69) 



(a 2 + Aa 2 + ^(C?- X D °C^ (70) 



The CP asymmetry can be obtained using the relation 

\A\ 2 - \A\ 2 _ 2\T SM+H \\E SM+H \sm{<)) 1 -0 2 )sin(-q g ) 

with <fii = Arg[af M+H ] and = Arg(xE)- As an example let us take Re(e^ 2 ) = 0.04, 
Im(e2 2 ) = 0.03 which is allowed point for tan/3 = 10. In this case we find that for a value 
of m H ± = 500 GeV we find that A CP ~ -3.7 x 10~ 5 while for m H = 300 GeV we find that 
Aqp — — 1 x 10~ 4 . Let us take another example where Re(e2 2 ) = —0.1, Im^e^) = —0.3 
which is allowed point for tan (3 = 500 and m H ± = 300 GeV. Repeating the same steps as 
above we find that Acp — 5.3 x 10~ 2 . Clearly in charged Higgs models the predicted CP 
asymmetry is so sensitive to the value of tan (3 and to the value of Higgs mass. 



V. CONCLUSION 

In this paper, we have studied the Cabibbo favored non-leptonic D° decays into K~n + . 
We have shown that the Standard Model prediction for the corresponding CP asymmetry is 
strongly suppressed and out of experimental range even taking into account the large strong 
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phases coming from the Final State Interactions. Then we explored new physics models 
taking into account three possible extensions namely, extra family, extra gauge bosons within 
Left-Right Grand Unification models and extra Higgs Fields. The fourth family model 
strongly improved SM prediction of the CP asymmetry but still the predicted CP asymmetry 
is far of the reach of LHCB or SuperB factory as SuperKEKB. The most promising models 
are no-manifest Left-Right extension of the SM where the LR mixing between the gauge 
bosons permits us to get a strong enhancement in the CP asymmetry. In such a model, it 
is possible to get CP asymmetry of order 10% which is within the range of LHCB and next 
generation of charm or B factory. The non-observation of such a huge CP asymmetry will 
strongly constrain the parameters of this model. In multi Higgs extensions of the SM, the 
2HDM type III is the most attractive as it permits to solve at the same time the puzzle 
coming from B — >■ tv and give a large contribution to this CP asymmetry depending on the 
charged Higgs masses and couplings. A maximal value of 5% can be reached with a Higgs 
mass of 300 GeV and large tan/3. 
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Appendix A: Operators and other definitions 

We start by defining X^p 3 , where Pi denotes a pseudoscalar meson, as follows 

x^ = ifpA 2 P2P .X 2P3 ( m \) ( A1 ) 
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where Ap 2?3 = mp 2 — mp g . In terms of Xp^ p we find that 



< n + \u lfll5 d\0 >< K-\s^c\D° > = -X£ 0K - 
< K~7r + \s^d\0 >< 0\u W c\D° > = 



(A2) 



Using Eq. ( 1A2I) we get 



< K-n + \O l \D Q > = < K-n + \s lllCL u lfl d L \D° >=< 7r + |u 7 ,A|0 X K-\s 1iiCl \D° > 

< ir-vr+IOsiD > = < X-7r + |M7 M c L S7AI^° >=< ^ + |s 7a A|0 X 0|m7 m c l |D° > 

+ i < 7r + |u 7/ A|0 X X-I^ID >= + ijCSJjf- 

< ^TT+I^C/j^rf^lD > = < TT+l^rf^lO X K-|s7 M Cjj|Z) > 

+ -jy < K-7r+\s^d R \0 X 0\u lt ,c R \D° >=-< K-n+ldlD > 

< K~7c + \u lfl c R s-f^d R \D° > = < K-jT+ls-f^dRlO >< 0\uj^c R \D° > 

+ ^ < vr + |M7^|0 >< K-\s^c R \D° >=-< K-7i+\0 2 \D > 

< if-vr+I^CL^G^lD > = < 7r + |M7 M rf fi |0 X K-\s ltl c L \D Q > 

- 1 < X-7r + |^ s+ p|0 X Oln^plD >= - ^X°°X^ + 

< K-7T + \u-f^c L s-f^d R \D° > = < K~7r + \sj fl d R \0 >< 0\uj^c L \D° > 

- ^ < 7T+|ud 5+ p|0 X ^"l^plD >= -Xf^ + + 

< ir-vr+l^cp^^lD > = < tt + |m7^ l |0 X K-\s^c R \D° > 

- |f < X-vr + |^ s _p|0 X OIwcs+pID >= X$ K - + 

< Jsr-vr+l^cp^rfLlL) > = < K~7r + \s^d L \0 >< Ol^cplD > 

- 1 < n + \ud s , P \0 >< K-\7sc s+P \D° >= X^ n+ - l x «* 'X& K - 

(A3) 
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and for the scalar ones 
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(A4) 

where the Fierz's ordering has been used 

#1^2, 5±P^3^4, S±P = -2^i^4, S±P^2, 5±P ~ ± 75 ) ( 1 ± 75)<V02 

1 — 1 1 

2(T a ) Q/3 (T a ) 7(5 = 5 a s5p y — — 5 a /3#7<5, (T a ) a/ 3(T a ) 7 5 = 5 a s5p y ~ (T a )as(T a )i3j 

(A5) 
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